Finite element method for solving neutron transport problem in two-dimensional cylindrical geometry

堀上 邦彦; 中原 康明; 藤村 統一郎; 大西 忠博*

JAERI-M 5793, 48 Pages, 1974/07

JAERI-M-5793.pdf:1.1MB

2次元(r、Z)体系での中性子輸送方程式を有限要素法をにより解くアルゴリズムを開発した。有限要素法は空間変数に対してのみ適用し、角度変数に対してはS法を用いた。(r.Z)平面を幾つかの長方形に分割し、それぞれの長方形の上でラグランジュ補間多項式を前もって作っておき、角度依存の中性子束をそれらの一次結合で表現する。一つの長方形の上で定義されるラグランジュ多項式の数は4、8、9の場合を考慮し、多項式の次数はr、Zの双一次、三次、双二次をそれぞれ対応させた。連続解を得るアルゴリズムと不連続解を得るアルゴリズムとを分けて説明し、両者いずれの場合においても、結合係数を決めるために適当な剰余を定義し、ガレルキン法を適用した。また連続解を得る解法の一つとして、中性子の保存則を表わす式を解くアルゴリズムも示した。